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Abstract 

High-precision interpolation of LISA phase measurements allows signal reconstruction and for- 
mulation of Time-Delay Interferometry (TDI) combinations to be conducted in post-processing. 
The reconstruction is based on phase measurements made at approximately 10 Hz, at regular in- 



^ I tervals independent of the TDI delay times. Interpolation introduces an error less than 1 x 10 

with continuous data segments as short as two seconds in duration. Potential simplifications in the 
^^ . design and operation of LISA are presented. 
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I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) is a mission to detect gravitational 
waves in the frequency band from 0.1 mHz to 1 Hz. The LISA constellation consists of three 
spacecraft flying in a heliocentric, Earth-trailing orbit, with separations oi L ^ 5 x 10^ m. 
Each spacecraft contains two proof masses that are shielded from external disturbances. 
To detect a passing gravitational wave, the change in separation 6L of the proof masses in 
different spacecraft must be monitored with a precision of SL/L < 10~^°/vIIz using laser 
interferometry. This fractional length stability is far better than the fractional frequency 
stability of the laser source, which is expected to be 5u/u > 10^^'^/vIIz. Degradation in 
sensitivity due to laser frequency noise could be avoided by operating the constellation as 
a Michelson interferometer with equal armlengths. Unfortunately, the orbital dynamics of 
the constellation make it impracticable to equalize the LISA armlengths accurately enough 
to cancel the excess frequency noise. Time delay interferometry (TDI) [l| is a technique to 
remove the otherwise overwhelming laser frequency fluctuations. TDI cancels laser frequency 
noise by combining phase measurements made at different times. The required timing of 
the measurements is set by the light travel times between the LISA spacecraft, and it must 
be accurate to 100 ns to meet the laser frequency noise suppression requirements [2|. 

One obvious method to achieve this timing accuracy is to measure the phase with a 
10 MHz sampling frequency. Selecting the nearest-neighbor samples would then provide 
the requisite 100 nsec timing resolution. This approach, however, would require data to be 
transmitted between spacecraft or back to Earth at the rate of approximately 10^ bits/s. 
The current design for TDI is to sample the phase at a much lower data rate, in the range of 
2 to 10 Hz, with 100 nsec accuracy triggering of the phasemeters J2|, yj . This approach poses 
a number of technical challenges. To ensure a timing accuracy of 100 ns, the absolute lengths 
of the arms must be known to an accuracy of 30 m when the measurement is made. For some 
TDI combinations, each spacecraft must have knowledge of its nonadjacent arm's length. 
Also, the clocks on different spacecraft must be synchronized at the 100 ns level. Errors in 
armlength knowledge or clock synchronization would lead to an irreversible corruption of 
the TDI combinations. 

An alternative approach is to sample the phase with a low rate at equally spaced times, 
and to reconstruct the phase at intermediate times by interpolation. Interpolation must 



be implemented with exceptional accuracy for effective cancellation of laser frequency noise 
by subsequent TDl processing. Tinto and colleagues j^ examined one possible method 
of interpolation and found that months of uninterrupted data around the time of interest 
are needed to achieve the necessary accuracy. This implies that months of data would be 
unusable at the beginning and end of a measurement, and levies extreme requirements on 
instrument reliability and operating duty cycle. The interpolation technique was deemed 
infeasible, and the triggered measurement approach was adopted. 

In this article we demonstrate that interpolation is feasible and that it can produce the 
required accuracy with less than two seconds of data. We discuss the significant simplifica- 
tion in the design and operation of the LISA mission resulting from this change. The method 
is based on fractional delay filtering j^ , a mature technique in digital signal processing. 

II. INTERPOLATION BY FRACTIONAL-DELAY FILTERING 

We specify that the interpolation error be less than 1 x 10~^ cycles/ a/Hz for frequency 
components from 1 mHz to 1 Hz. This noise level is approximately a factor of 10 below the 
phase noise contribution of shot noise. Below 1 mHz the requirement is relaxed, as the 1//^ 
proof mass displacement noise dominates shot noise and a larger interpolation error can be 
tolerated. Assuming that the laser frequency noise produces approximately 100 cycles/-\/Hz 
at the phasemeter output, interpolation must have a fractional error of less than 1 x 10~® 
for frequency components in the 1 mHz to 1 Hz range. 

We assume that the LISA phase measurements will be recorded with a /^ = 10 Hz sam- 
pling rate. The sampling rate must be high enough to accurately reproduce the phase infor- 
mation in the LISA signal band, and to avoid adding noise from aliasing of higher-frequency 
phase noise, at the 10~^ cycles/ vHz level. Moreover, the performance of the interpolation 
schemes considered below improves with oversampling. Ultimately, the sampling rate will 
be determined by filtering requirements on the phasemeter and by the availability and cost 
of telemetry bandwidth to Earth. 



A. Perfect interpolation and fractional-delay filters 

Interpolation is tire process^ot reconstructing the amplitude of a regularly sampled signal 

between samples. Shannon [5] proved that a bandlimited signal sampled at a sufficiently 
high frequency can be reconstructed perfectly by convolving the discrete time series with a 
continuous sine cardinal function sinc(/st) = sin^Ti fgt)/ {it fst). 

Sine interpolation can also be viewed as applying an acausal finite-impulse-response (FIR) 
filter to the sampled time series. The filter kernel (impulse response) is a sampled version 
of the sine function. In effect, instead of interpolating the signal we interpolate the filter 
kernel. As the sine function is a known analytic function, the filter kernel can be interpolated 
with arbitrary accuracy simply by time shifting the argument of the sine. In general, the 
interpolated signal s{n — D) is the discrete convolution of the original signal s{n) with the 
shifted kernel: 

s{n-D) = s{n)*h{n-D), (1) 

where n is the sampling index, D is the delay in samples (— | < -D < |) and h{n) is the 
filter kernel. For sine interpolation, h{n — D) = sinc(n — D). 

With zero delay (D = 0), sine interpolation corresponds to a FIR filter with delta function 
impulse response [see Fig. ^a)], since for integer n smc{n)=6no where 6nk is the Kronecker 
delta function. li D y^ 0, we obtain a FIR filter with a non-delta impulse response [see 
Fig. E^b)], which has the effect of applying the fractional delay D to the original time 
series. Errors in fractional-delay filtering are caused by the finite length approximation of 
the infinitely long delayed-sinc filter. 

B. Truncated-sine fractional-delay filters 

The simplest finite-length approximation to the ideal delayed-sinc filter is obtained by 
truncating the kernel. Filtering by a truncated-sine of kernel length N can be written as 
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D)= Y^ s{n + k) smc{D - k) (for odd A^). (2) 
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In the following discussion we restrict ourselves to filters where A^ is odd for simplicity. 
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ipugh for a given filter order A^ the truncated-sine is optimal in a least-squares 
, its frequency response is far from ideal (unity magnitude), exhibiting significant 
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FIG. 1: Sine FIR filter kernel values (circles) with delay set to (a) D = 0, and (b) D = 0.3. 



ripple even at low frequencies. This is unacceptable for TDI, where very high fidelity is re- 
quired in the 1 mHz to 1 Hz measurement band. In fact, Ref. J2| showed that truncated-sine 
interpolation becomes sufficiently accurate only for very large N. Figure |21 shows the inter- 
polation error versus A^ for several filters, including the truncated-sine. The interpolation 
error e is defined as the maximum difference of the filter's frequency response and the ideal 
frequency response, q~^'^'^S^iS^ for frequencies between 1 mHz and 1 Hz: 



e = max(|iJ(/) — e 
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where H{f) is the Fourier transform of h{n). We used D = 0.5, which is expected from 
theory to be the worst case. With truncated-sine interpolation, e ~ 1/A^. Sampling at 10 Hz, 
this filter would require a kernel almost four months long, A^ > 10^, to achieve e < 10^*^. 
This means that two months of data at the beginning and end of each measurement period 
would be unusable. 




Kernel length A'^ 

FIG. 2: Comparison of interpolation error for four interpolation methods; truncated-sine filter 
kernel, Blackman windowed-sine filter kernel, Blackman^ windowed-sine filter kernel, and Lagrange 
interpolation. 

C. Windowed-sine fraetional-delay filters 

The ripple in the frequency response of the truncated-sine filter can be significantly 
reduced by windowing the filter. 
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D) = y^ s{n + k)w{k)smc{D — k), 



(4) 



where the window w{k) goes smoothly to zero for k = ±{N — l)/2, so that the endpoints are 
tapered to zero instead of abruptly truncated. Of several conventional windows jGj] tested, 
the Blackman function 
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was best, producing e < 10"^ for A^ > 345 (see Fig. 12)) corresponding to a loss of 14.4 seconds 
of data at the beginning and end of each measurement period. In comparison, the TDI 



combinations need several L/c arm travel times, or at least 65 seconds, to gather enough 
data to cancel laser noise. As seen in Fig. £ ~ 1/A^^ for Blackman windowed- sine filters. 
One simple modification to the Blackman windowed-sine filter kernel is to apply the 
Blackman function more than once. Our tests showed that using wl{n) (applying the Black- 
man three times) produced e < 10~* for A^ > 21, corresponding to a loss of 1.1 seconds of 
data at the beginning and end of each measurement period. 

D. Lagrange filter 

A more accurate filter at low frequencies can be found by requiring a maximally flat 
frequency response at DC |^. This filter kernel is equal to the Lagrange polynomial [3] 
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"-(") = n ^. (6) 



_ JV-1 
2 



fc- 2 



where tn = ^^^^ + D. The Lagrange filter can be compared to filters in Section III CI by 
expressing its kernel as a windowed-sine function (Eq. |3)), with the window 

"^^"^ = ^ (!n) (n + ^"-1)72) ' (^) 

where the binomial coefficient is extended to non-integer arguments by the generalized fac- 

n 

torial function (Gamma function) ^]. Figure El shows the Wb{n), wf{n) and wiin) with 
A^ = 31. 

The performance of the Lagrange filter for this application is excellent, meeting the 
requirements with A^ = 15 (1.5 seconds) as shown in Fig. |21 Although Lagrange interpolation 
is known to produce large spurious oscillations at the ends of the interpolation interval, this 
does not occur when the kernel is fully immersed in the signal. By accepting only data 
where the kernel is completely immersed, we obtain excellent performance at the expense of 
losing A^/2 points at the beginning and end of each measurement period. 

Lagrange interpolation is related also to the Thiran infinite-impulse-response (IIR) 
fractional-delay filter ^, which has nominally flat frequency response. The performance 
of the Thiran filter in our application is comparable to the Lagrange windowed-sine FIR 
of the same order; but the latter is favored on grounds of simplicity, especially for time- 
dependent delays. 
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FIG. 3: Blackman, Blackman^ and Lagrange windows with A^ = 31. 

The filters were tested both by interpolating known analytic functions and interpolating 
bandlimited white noise. The bandlimited noise was generated with a sampling rate of 
10 MHz and resampled at times t = n/ fs, fs = 10 Hz. The 10 Hz signal was delayed by D 
samples and compared to the original 10 MHz signal resampled at times t — D / fg. The test 
results agreed with the calculated fractional error shown in Fig. |21 

E. Further tests 

We have characterized the error of fractional-delay filtering with fixed delays. In orbit, 
the delays will slowly vary due to the changing arm lengths, and so we also tested Lagrange 
filtering with varying delay. For this test, we generated a time series of white noise, ban- 
dlimited to 2.5 Hz and sampled at 10 Hz (oversampling factor of 2); we then used Lagrange 
filters of increasing order to interpolate the noise to the original sampling times shifted 
by delays ranging linearly in time from D = —0.5 to D = (no delay), during a period 
of 5 X 10^ seconds. This arrangement approximately simulates the slow variation in the 
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LISA armlengths (which determine the TDI-mandated delays). Figure 0] shows spectra of 
the interpolation error, along with the spectrum of the original white noise. The required 
interpolation accuracy is achieved at all frequencies in the measurement band for A^ > 16 
(window length of 1.6 s). 

If the requirements become more stringent, for example due to increased laser frequency 
noise, it is easy to improve the performance to the desired level by increasing the length of 
the filter kernel. By setting the filter to A^ = 31 (3 seconds in length at 10 Hz), e can be 
reduced to 1 x 10"^^. 



10^ 



N 

I 10-2 
o 

o 10-^ 

iij 

c 
o 

8. 10"' 

B 
c 



1 cycle/rt(Hz) noise band limited to 2.5 Hz 



0.4 s window 



10"' 



ft. 8^ window" 



1.2 s window 



10 



10 



1.6 s window 



2.0 s window 



WHMMf"'" 



end of measurement 
band 




10' 



10 
Frequency [Hz' 



FIG. 4: Amplitude spectral density of interpolation error for Lagrange filters, shown with the 
spectral density of the initial 2.5 Hz-bandlimited noise. Spectral density is estimated by a triangle- 
windowed, averaged periodogram. 
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FIG. 5: Real-time TDI, in which the delays must be known at the time of measurement, and 
the TDI combinations are computed before transponding to Earth. TDI combination X{t) shown 
(Eq.lHl). 

III. IMPLICATIONS FOR LISA 



To illustrate the design simplifications enabled by high-precision interpolation, consider 
the simplest TDI combination-the Michelson combination X(t) jjj: 



X{t) = [S2i(t) - S3l(t)] - [S2i(t - 2L3/C) - S3l(t - 2L2/C)]. 



(8) 



where Smi is the phase measurement made at Spacecraft 1 of the light received from Space- 
craft m, Ln is the length of the arm opposite Spacecraft n, and we are assuming that the 
six LISA lasers are phase locked j^ . The first two terms represent the optical phase of two 
arms of a Michelson interferometer, and the second two terms represent the same quantity 
with the specified delays. 

The implementation of TDI based on timed triggering J2], which we designate "real- 
time TDI," calls for all four terms in Eq. (jHI) to be explicitly measured, combined on- 
board, and sent to the Earth. This is illustrated in Fig. for the measurement of X{t) 
aboard Spacecraft 1. The blocks labeled pd2i and pd^i represent the radio-frequency beat 
signals from the photodiode outputs, containing laser frequency noise superimposed on the 
gravitational wave signal. The Ranging block represents the system that measures distances 

10 



Pd2i 



l*> PM I- 



S2l(t) 



pdn 



I*- PM I ►" 



Earth 



■ _ Interpolation C TDI ■— ► 



S3l(t) 



FIG. 6: Post-processed implementation of TDI combination X{t). 

between spacecraft, and computes the delays 2L^/c and IL^jc required to assemble X(t). 
The telemetry inputs to the ranging blocks (not shown) contain ranging data measured on 
Spacecraft 1 and Spacecraft 2. Phase measurements are made by the phasemeter (PM) 
blocks. In Fig. the phase signals are delayed electronically by the Delay blocks, which 
implement a variable delay as controlled by the ranging system. Equivalently, identical 
phase signals can be fed to two PM blocks, and the timing of the phase measurements 
can be set by adjustable triggers from the ranging system outputs. The final output X{C) 
is free of laser frequency noise for fixed armlengths. For time-dependent armlengths, we 
expect J9| that velocity-correcting or "second generation" TDI combinations will be required. 
They have roughly double the measurements of length-correcting or "first generation" TDI 
combinations such as that shown in Fig. El 

As we have demonstrated in this paper, the delayed phase measurements required for 
TDI can alternatively be inferred by interpolation from an equally-spaced sequence at a 
relatively low rate. This implementation is referred to as post-processed TDI and is shown 
in Figure El The elimination of the need for ranging knowledge at the time of measurement 
simplifies the implementation, as is evident comparing Figs. JHl and JHl Ranging information 
will still be needed as input to signal reconstruction, but it can be transmitted to Earth 
independently of the phasemeter signals. Alternatively, post-processed TDI can be imple- 
mented without explicit ranging by determining the delays using autocorrelations [3], or 
by adjusting the delays in post-processing for minimum sensitivity to laser frequency vari- 
ations. With post-processed TDI it is no longer necessary to synchronize the clocks on 
different spacecraft. Clock synchronization error can be corrected simply by time-shifting 
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the data in post-processing. The decouphng of ranging from phase measurements and the 
removal of the need for clock synchronization allows reduction-or possibly elimination-of 
inter-spacecraft communications. 

Post-processed TDI also allows complete flexibility in combining phasemeter signals. All 
the raw data are available for processing by any TDI algorithm, including ones not developed 
until after the data are in hand. The delays can be adjusted to optimize the suppression of 
laser frequency noise; by contrast, if there is an error in triggering in real-time TDI, noise is 
irrevocably added. Post-processed TDI simplifies the phase measurement hardware, allowing 
all possible TDI combinations to be constructed from one constant-rate phase measurement 
per photodetector. 

A significant operating cost of LISA will be telemetry to Earth of science data. Real-time 
TDI requires one data stream per TDI combination. Post-processed TDI requires one data 
stream per phasemeter; more precisely, per phasemeter that does not have its output held 
fixed by a high-gain control system. We expect that post-processed TDI will require fewer 
telemetry signals than real-time TDI, but this depends on details of hardware design and 
on data requirements that are currently under consideration. 

Other factors influencing overall telemetry costs are the data rate and number of bits per 
datum. The ultimate LISA data will have a signal bandwidth of 1 Hz and a dynamic range 
large enough to encompass both the largest expected gravitational wave signal (or perhaps 
the largest instrumental effect) and shot-noise limited sensitivity. Each output of the real- 
time TDI signal chain {X(t) for example) will have essentially eliminated laser frequency 
noise before data are transponded to Earth, reducing the dynamic range requirement. A 
nominal datum size is 20 bits per sample. The data rate for each TDI combination is 
2 samples/second, in keeping with the 1 Hz requirement and the Nyquist limit. 

In comparison, post-processed TDI must transpond large laser frequency noise superim- 
posed on the small gravitational wave signal. Larger dynamic range is required, perhaps 
30 bits per sample. The sample rate for post-processed TDI is likely to be greater than 
2 samples/second for two reasons. Firstly a larger sampling frequency may be needed to 
provide the more stringent anti-aliasing filtering needed when laser frequency noise is present. 
Secondly an oversampling factor may be needed for the interpolation procedure. The min- 
imum sample rate imposed by post-processed TDI is still under study; our initial estimate 
is 10 samples/second (oversampling factor of 5) or less. 
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Although the interpolation algorithms presented here may not be optimal, they demon- 
strate the feasibility of post-processed TDI. They serve as a proof of principle, and provide 
guidance for the design of LISA with significant simplification in several respects over real- 
time TDI. 
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